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Abstract 

In a series of publications [1,2], Liischer et al. have demonstrated the usefulness of the 
Schrodinger functional in pure SU(2) and SU(3) gauge theory. In this paper, it is shown 
how their formalism can be extended to include fermions. In the framework of Wilson's 
lattice QCD, we define the Schrodinger functional by making use of the transfer matrix 
formalism. Boundary conditions for the fermions arise naturally. We then take the naive 
continuum limit of the action and show that no lattice peculiarities are left over. The 
corresponding free Dirac operator has a unique self-adjoint extension with purely discrete 
spectrum and no zero modes. 
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1. Introduction 



Recently, it has been shown how the Schrodinger functional can be used for a the- 
oretical determination of the running coupling in asymptotically free gauge theories [1]. 
The subsequent applications to the SU(2) and SU(3) Yang-Mills theories have been very 
successful [2], so that one may feel encouraged to extend this framework to QCD. 

The Schrodinger functional as defined in reference [1] is the euclidean propagation 
kernel for going from an initial field configuration at time xq = to a final configuration 
at euclidean time xq = T. Following Feynman, such a transition amplitude can be written 
as a path integral with the usual (euclidean) action, where the integral runs over all 
interpolating field configurations. The quantum fields then fluctuate around the "classical 
path", i.e. a classical background field, which may be modified by tuning the initial and 
final field configurations. 

The Schrodinger functional is important for the application mentioned above, since 
it allows to define a running coupling as the system's response to the induced background 
field. This definition is not inherently based on the perturbation expansion and is therefore 
also well-suited for computations in the non-perturbative low-energy domain. 

It is the purpose of the present paper to show how the Schrodinger functional can be 
defined in QCD. In a theory with fermions, it is not clear what one should take as initial 
and final field configurations. In the euclidean framework this amounts to the problem of 
finding adequate boundary conditions for the quark fields. The difficulty arises because 
the Dirac equation is only of first order, and therefore, one can only prescribe the values 
of half of the spinor components at the boundary. 

To clarify the situation, it is most convenient to define the Schrodinger functional 
on the lattice, using the transfer matrix formalism. The lattice action and the boundary 
conditions then emerge naturally. As a first check, we then take the continuum limit of the 
lattice action and verify whether this leads to an acceptable classical continuum theory. 

The paper is organized as follows. Section 2 introduces some of our conventions. The 
case of pure gauge theory is reviewed in section 3, stressing the role of the transfer ma- 
trix. We then introduce the transfer matrix for Wilson fermions and discuss a convenient 
technique for obtaining its Grassmann representation (sect. 4). The Schrodinger functional 
is defined and the lattice action and the boundary conditions are read off from its path 
integral representation. Sects. 5 and 6 are concerned with the classical theory, and we end 
with a few remarks concerning questions of immediate interest. 
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2. Basic conventions 



In order to introduce some conventions and notations, we consider QCD on a finite 
hyper cubic euclidean lattice 

x/aeTZ*, < x° < T, < x k < L, k = 1, 2, 3} . (2.1) 

A classical quark field is an assignment of a Dirac spinor ^(x) with components ^(x) saa 
to each lattice site. The indices s and a denote the flavour and colour degrees of freedom 
and a is the Dirac index. A gauge field is an ensemble of matrices U(x,fi) £ SU(iV), 
residing on the bonds (a;, x + afi) of the lattice (fi denotes a unit vector in /x-direction). 
Throughout the paper, periodic boundary conditions in the spatial directions are assumed 
for all fields. Assuming periodicity in the time direction too, one may write down the 
action, S = S g + Sf, according to Wilson [3] as follows 

S M=\ E E trjl-P^)}, (2.2) 

3 

Sfbl>, j>, u] = a 3 E - E $( X M X > aO£(i - iM x + «A) 

+ ${x + afi)U(x, a*) _1 /C(1 + lM( x )] }• ( 2 - 3 ) 

Here, P fJil ,(x) = U(x,fi)U(x + afi,v)U{x + av, fi)~ 1 U(x, v)~ l is the usual plaquette field 
and K. is the hopping parameter matrix, assumed to be diagonal and non-negative, its 
relation to the bare quark mass matrix m being K~ x = 2ma + 8. Our 7-matrices are 
hermitian and satisfy 

{7^,7^ = 2^. (2.4) 

In the continuum one is used to a different normalization of the spinor fields, namely 
ipc = (2/C) 1 / 2 ^ (c stands for continuum). Later on, we will consider the continuum theory 
and then use this normalization. If one introduces the covariant lattice derivatives 

V^(x) = -[U (x, n)if>(x + afi) - if>(x)] , 

\ (2.5) 
V^z) = - [${ x ) ~U(x- afi, fi) 1 tp(x - afi)] , 

and the difference operators V M = (V^ + V*)/2 and A = — X^=o ^aj^^i one can rewr ite 
the fermionic action Sf in the form 

S f = a i E &(s){E^' + ro +^ aA }^0 B )> (2-6) 

xev E a 4 

which is reminiscent of the euclidean action in the continuum. 
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3. Pure SU(iV) gauge theory reviewed 



The Schrodinger functional in pure SU(iV) gauge theory has been discussed in refer- 
ence [1], where also a formal approach in the continuum can be found. We find it useful to 
go through the lattice definition of this functional in some detail. First of all, it is evident 
that the pure gauge theory functional will constitute an integral part of the full QCD 
functional. Furthermore, from the discussion of the pure gauge theory, it will become 
clear how to proceed for lattice QCD. 

The starting point is the quantum mechanics of the pure gauge theory. One considers 
^-periodic gauge fields on a spatial lattice T, which can be obtained as a sublattice of the 
euclidean lattice Te (2.1) for a fixed value of xq (xq = for instance). Wave functions are 
defined as the complex-valued square integrable functions of the gauge fields. They form 
a Hilbert space with the scalar product given by 

<^ 2 >= / V[U]<pl(U)<p 2 (U), V[U]= J] JJdtf (*,*)■ (3-1) 

J xerfe=i 

The wave function <p is the Schrodinger representation of a state vector \<p). It can be 
obtained from the latter by introducing a complete set of states \U), labelled by classical 
gauge fields U such that 

(U\<p) = <p(U) . (3.2) 

Since the field operators £7(x, k) are diagonal in the Schrodinger representation, the states 
\U) are required to be eigenstates of £7(x, k) with the corresponding classical field as 
eigenvalue. Completeness implies the decomposition of the identity operator 

1 = J V[U]\U)(U\, (3.3) 

which allows to retrieve a state vector from its Schrodinger representation. 

Only gauge invariant wave functions, which satisfy f{U) = <p(U v ) are physical. Here, 
U v denotes the gauge transformed field U v (n, k) = F(x)£7(x, fc)F(x + afc) -1 , with F(x) 
being an ^-periodic SU(iV)-matrix field on T. Any wave function can be projected on the 
physical subspace with the help of the projector IP, 

P<P(U)= I Y[dV(x)<p(U v ). (3.4) 
J xer 

The (euclidean) time evolution of a wave function is now provided by the transfer matrix 
T = To IP, which can be thought of as the euclidean step evolution operator e _IHa (with 
Hamiltonian IH). The operator To acts on wave functions as an integral operator 

T <p{U) = (U\T \<p) = J V[U'] (U\T \U')<p(U'), (3.5) 
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with kernel K [U, U'\ = (Z7|T |Z7') = exp -AS[U, 1, U'\. When making contact with the 
path integral formulation, one usually imagines a double layer of two space lattices T, 
separated by a euclidean time a. Apart from the factors 1/2, the function AS[U, V, U'] = 
S p [U]/2 + Sk[U, V, U'] + S p [U']/2 then equals the euclidean action of this double layer, 
provided V is interpreted as the temporal component of the euclidean gauge field. For the 
Wilson plaquette action (cf. sect.l), one defines 



^[P] = iE E tr{l-P fcI (*)}> 

Grfe 3 =1 (3-6) 

S k [U, V, U'] = \ Y, E tr < 2 - Pfe °( x ) - ^(x)" 1 } . 
9o 



xer fe=i 



For obvious reasons, the operator To is referred to as transfer operator in the temporal 
gauge. Using the gauge invariance of the action, 

AS[U, V, U'] = AS[U, 1, U^' 1 ] = AS[U V , 1, U'], (3.7) 

and the invariance of the measure, one easily sees that To commutes with IP and the 
integral kernel of the transfer matrix T = To IP is given by 

(U\T\U') = / II dF ( x ) e" A5[C/ ' F ' U '^. (3.8) 
J xer 

The role of the projector IP therefore is to supply the temporal components of the euclidean 
gauge field. 

Euclidean time transport over a time T = Na is now obtained by iterating the action 
of the transfer operator N times, viz. 

{T N <p){U N ) = J V[U ] (U N \T N \U )<p(U ). (3.9) 

The kernel can be evaluated by inserting the identity (3.3) N — 1 times, 

. N-l 

(U N \T N \U ) = J J] P[^](C/ J v|T JV |C/ J v-i)---(C/ 2 |T|C/ 1 )(C/ 1 |T|C/o). (3.10) 



i=i 



If one inserts the expression (3.8) and identifies the index I with the euclidean time x^/a 
one obtains the path integral representation of the integral kernel 



'i-iriv 



(W'\T 



\W) = j B[U]e- S \- U \ (3.11) 
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where the integration is over the euclidean gauge fields U(x, /x) with < xq < T for /x = 
and < xq < T for /x = 1, 2, 3. The action S, given by 

s w 4EE E™( p ) tr ii - p»a*)}, (3.i2) 

ffo xer x °=o 

is the usual Wilson plaquette action if it were not for the weight factor w(P). This factor 
is set to 1 for all plaquettes P except the space-like ones at the boundary (x$ = , T) 
where it is set to 1/2. Furthermore, the euclidean gauge field is subject to the boundary 
conditions 

U(x,k)\ x o =0 = W(x,k), and U(x, k)\ x o =T = W'(x, k) . (3.13) 

It is the integral kernel (W|T T / a |W A ) which is referred to as Schrodinger functional. How- 
ever, being interested in the continuum limit, we will again follow reference [1] and con- 
sider the Schrodinger functional as functional of the continuum gauge potentials Cfe(x) 
and C£(x). The relation with the lattice boundary fields W and W is then established 
through parallel transport along the lattice bonds, i.e. 



W(x, k) = Vexp^aj^ dtC k (x + ak - tak) j, (3.14) 

and analogously for the primed fields (the symbol V means a path ordered exponential 
such that the fields at larger values of t come first). The Schrodinger functional is then 
written as follows 

Z[C, C] d = (W'\T T/a \W) . (3.15) 



4. Extension to QCD 



In order to extend the definition (3.15) to QCD, one may use the same strategy as 
in the pure gauge theory case. This means, one starts with the transfer matrix, raises 
it to the iV-th power and represents it as an integral operator, acting on suitable wave 
functions. The integral kernel of this operator will then be the Schrodinger functional 
one is looking for. The transfer matrix for lattice QCD with Wilson fermions has been 
constructed by Liischer [4] some time ago, when he proved positivity of Wilson's lattice 
QCD. Except for some details and different notations, we will follow this approach. 
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4-1 Transfer matrix 

The starting point is again the time zero Hilbert space which is now the product space of 
the pure gauge theory Hilbert space and a fermionic Fock space, generated by the action of 
a set of canonical operators %(x) saa , X^(x) saa on the vacuum state |0) . These operators 
reside on the sites of a spatial lattice T with periodic boundary conditions. The flavour, 
colour and Dirac indices are the same as for the classical fields (cf. sect. 2). By definition, 
one has 

{x( x ) S a«> tfWtbp} = SstS ab 5 af 35(x - y) . (4.1) 

In fact, it is an important observation that one should not impose canonical anti- 
commutation relations on the operators if) and ^ which correspond to the fields appearing 
in the classical action. Instead, these are related to the canonical operators through 

(4.2) 

^( x ) M «= X) * f (yW 5 1/2 (y' x W 

y,t,b,P 

and thus fulfill the anti-commutation relation 

{^W.aa»^ f (yW} = ^B-\^y) saa , m . (4.3) 

The matrix B is at this stage only required to be hermitian and positive. When the 
connection with the euclidean formulation is established, one finds f 

3 

5 ( x > y)s aa ; tbp = <W { S tt 8 ab 8(x, y) - A^V [U (x, k) ab 8(x + ah, y) 

1 t^i (4.4) 

+ U(y,k); b 1 5(y + ak,x)}}. 

The matrix B can be obtained as the kernel of the difference operator 5 = 1 — 6K. — 
a 2 JC Y^k=i ^fe^fe' taking into account the periodic boundary conditions on the spatial 
lattice. It is a hermitian and positive matrix if all hopping parameters are smaller than 
1/6, which will be assumed in the following. 

A general state \X(U)) of the Hilbert space evolves in euclidean time through the 
action of the transfer operator 

T |X(10> = J V[U']T [U, U'] \X(U')) . (4.5) 



■j" Concerning the spatial 5-functions, we distinguish between 6(x,y) with 5(x,x) = 1, and the 3- 
dimensional lattice 5-function 5(x — y) = a -3 6(x,y). 
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With respect to the gauge field Tq is an integral operator whose kernel has the structure 



T [U, U'] = f F (U)K [U, U']f F (U') . (4.6) 

The kernel Kq[U,U'] is the same as for the pure gauge theory [cf. eq.(3.5)], while the 
operator Tp(U) is given by 

f F (U) = det(2X:S) 1 / 4 exp(-*|(l - lo )C X ) exp(-^Mx) . (4.7) 

The adjoint spinor is defined as usual by % = jtflo an d the matrix M is related to B 
through M = 1/2 ln^^ZC) -1 ]. The anti-hermit ian matrix C is given by 



3 1 

C( X > y) saa . tb(i = S st Y (Tk)a/3 9 i U ( X > k )ab S (* + ak , y) 



k=l 



2 L v ' v ■ (4.8) 



U(y,k)-£6(y + ak,x)], 



and is obtained as matrix kernel of C = a Y^k=i Ik^k- 

In order to define the Schrodinger functional, one needs a representation of the transfer 
matrix as an integral operator acting on wave functions. It is well-known that any operator 
in Fock space may be represented as an integral operator in an associated Grassmann 
algebra. There is a nice technique to represent this isomorphism, based on so-called 
"Grassmann coherent states" (see e.g. [5]). To apply this technique, it is convenient to 
pass to a two-spinor representation and to choose a representation for the j- matrices. 
This is only an intermediate step and at the end one may re-express everything in a way 
that does not depend on the particular representation, chosen as follows (rj., k = 1, 2, 3 
denote the Pauli matrices) 

We will frequently use the projectors P± = |(1 ±70) which are diagonal in this represen- 
tation and used to project on the two-spinor operators, defined by 



X 



c ?t), X f = (^C). (4.10) 



The matrices B and C have projections B and C 

In the two-spinor space we introduce barred greek indices a, /3.., which take the values 1 
or 2. Classical spinor fields on the spatial lattice may be treated in parallel to the operator 
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fields, i.e. the classical fields x(x) and x(x) with 



X (x) = (i^VHx), x(x) = (^/ 2 )(x), 
are decomposed as follows 



x=( c +). x = (v + ,-0 



(4.12) 



(4.13) 



Whenever the four-component spinors \ an d X appear together with the two-component 
spinors tj, tj + and £, £ + within the same formula, this decomposition is tacitly assumed. 

4-2 Grassmann coherent states 

We proceed by introducing the states Grassmann coherent states (f/ + C + | an d \tjQ- These 
are Fock-states with Grassmann coefficients, such that the following holds 



(v + C + \fi\x) saSl = v + (x) saSl (v + C + \ 



(4.14) 



for all possible values of the indices (analogous relations hold for the operators £ and 
It is always assumed that canonical operators and Grassmann variables commute with 
each other. The Grassmann coherent states are eigenstates of the canonical operators, 
with the corresponding classical fields as eigenvalues. It is not difficult to see that they are 
uniquely determined if the normalization conditions (f/ + C + |0) = 1 = (0|r/C) are imposed 
(assuming (0|0) = 1). One then finds 



00 , u J(J-l)/2 



x,s,a,a 



1°)' (4.15) 



and, similarly, 



00 , ,J(J-l)/2 

(. + c + i=(oiE izlL ^ — 



j=o 



a3 X) {^( X ) S aa 7 / + ( X ) S a a + C(x) saa C + (x) saa } 



(4.16) 



x,s,a,a 



The Grassmann fields will become integration variables. When several integrations have 
to be carried out, we will need different sets of Grassmann variables which are associated 
to the same canonical operators. In this case we will use a prime or an index I with the 
understanding that e.g. 7/ ; (x) saC( will be the eigenvalue of 7/(x) saa when acting on the 
eigenstate |t/jC/)- I n the euclidean framework this label will then become the euclidean 
time index. 
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The isomorphism alluded to above takes a general Fock state \X) to its Grassmann 
equivalent, given by (t] + £ + \X). This is the wave function we were looking for. It only 
depends on the crossed Grassmann variables and therefore, such wave functions are some- 
times called "holomorphic states" in the literature [6]. 

Using the rules of Grassmann integration [7], one obtains the following decomposition 
of the identity operation in the Fock space 

w= /n n d ^ + ( x ) saa ^(x) saa dc + (x) saa dc(x) saa 

x£T s,a,a 

x ex P [-a 3 2 Ei^ + (y)^^(y)^ + c + (y)^c(y)^}] i^cx^+c+i - 

y t,b,p 



(4.17) 

It is this "completeness relation" which allows to recover a Fock state \X) from its wave 
function (t] + £ + \X), thus proving that one really has to do with an isomorphism. As can 
be seen from eqs.(4.1) and (4.14), both, fermionic operators and the Grassmann fields 
have canonical dimension 3/2 and the normalization of the Grassmann integrals is chosen 
as follows 



/ 



II II d ^ + ( x ) saa dr/(x) saa dC + (x) saa dC(x) saa 

xGrs ' a ' a n (4.18) 

exp [-a 3 {v + (y W v(y)tbf3 + C + (y)tbf3 C(yW} 



y t,b,p 



We are now in the position to express any Fock space operator f as integral operator in 
the Grassmann algebra by applying the isomorphism to the state T\X) and inserting the 
identity once 

(rj + C + \f\X) = J V[ X ']V[ X ']e-x'x'(ri + C + \T\Ti'C')W + C' + \X). (4.19) 

The measure appearing here is the same as in eq.(4.17) (cf. the remark after eq.(4.13)) 
and defined through 

V[X]V[ X ] = II II d x( x ) saa dx(x) saa . (4.20) 

x£T s,a,a 

Products and powers of operators can be treated similarly by inserting the identity oper- 
ation (4.17) between any two factors. 
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4-3 Grassmann representation of the transfer matrix 



All one has to do now, is to apply this procedure to the transfer matrix (4.6) and its iV-th 
power. Starting from equation (4.7), one first expresses the operator T F (U)Tp(U') in the 
two-component spinors and reorders it such as to have the crossed operators to the left of 
the uncrossed operators, 

fl(U)f F (U') = det(2A: J B) 1 / 4 det(2A: J B') 1 / 4 

X exp(^e- M Ce- M (i) exp(-^Mr) + (M^) (4.21) 
X exp(-^ t M'fj + CM'C f ) exp(Ce - *' C'e"*' rj). 

The corresponding integral kernel with respect to the Grassmann variables is then found 
to be 

( V + C + \fl(U)f F (U')\r)C) = det(BB') 1 / 2 exp(r ] + B- 1 / 2 2KCB- 1 / 2 C + ) 

X exp( V + B- 1 / 2 21CB'- 1 / 2 V - ^B'- X I 2 2KB~ X I 2 (+) (4.22) 
xexp((B'- 1 / 2 2KC'B'- 1 / 2 V ). 

By definition the crossed Grassmann variables tj + and £ + correspond to %P+ = x+ an d 
P-X = X- [ c f- eq.(4.13)]. Since the matrix B commutes with the projectors P±, one may 
also define a wave function X as a function of the components ^_ and viz. 

X(j + ,^,U) = (rj + C + \X(U)). (4.23) 

Performing the change of variables in the Grassmann integral, using eq.(4.22) and the 
transformation of the measure (4.20), ^[xjl^x] = T>\t\)\D\t\)\l det B, one finds the full 
action of the operator To on a wave function 

(T X)(i> + ,i{,_,U) = j V[U']V[iP']V[iP'] /det B' e-^' B '^' det(BB') 1 / 2 

x e -2{^P + KC^ - ^P+Kt/f' - $'P-Ktl> + ^'P-KC'^'} (4.24) 

X e-^^' 1 '^']^;,^,^')- 

As for the pure gauge theory, one introduces the projector IP on the physical subspace of 
gauge invariant wave functions, 

PX$+,tl>-,U)= I X[dV(x)X(i, + V-\Vil,-,U v ). (4.25) 
J xer 
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The operator To is again referred to as transfer matrix in the temporal gauge. The full 
integral kernel of the transfer matrix T = To IP = IP To in terms of the classical four- 
component spinors is now easily found, making use of the gauge invariance of the measure 

( V +( + \T[U,U']\ V ' (') = det(BB') 1 / 2 / II dF ( x ) e" A5[C/ ' F ' U ' ] 

*er (4.26) 



X e 



4-4 Definition of the Schr '6 'dinger functional 

It is not difficult to iterate the procedure and to evaluate the iV-th power of the transfer 
matrix. Its kernel 

(vM\V N [U N ,U ]\ Vo C )= / n V[UiMXiMXi]e~^ 

J 1=1 

(4.27) 

x (Vn(n\ t [ u n, ^iv-i]|f/jv-iCjv-i) • • • 

• • • (rf(t\nU2, U^n^in+CtinUuUo]^) 

can be obtained in terms of the classical quark fields t\) and tj) by inserting the expression 
(4.26) in the r.h.s. of this equation and changing the measure accordingly. One then 
obtains the following expression for the kernel of T^ 

f f N-l 

det(B N B o y/ 2 J Y[dV (x) J J] vmVtyWfa] JJdVi(x)e-^^ 



1=1 

N-l 



(4.28) 



x exp - J] {AS[U k+1 , V k , U k ] + ^ k P-2K C k ^ k - ^ k P-2K V k ^ k+1 

k=0 

- ^ k+l P + 2KV k ~ l ^ k + ^ fc+ iP + 2/CC fc+ i^ fc+ i} . 

As announced above, the label I is interpreted as euclidean time, i.e. one sets 

l = x°/a, N = T/a. (4.29) 

The link variables t//(x, k) and the fermi fields ^(x), ^(x) are then identified accordingly 
with U(x, k), ^(x) and tj>{x). The fields Vj, stemming from the projectors IP play the role 
of the time-like link variables V^(x) = U{x, 0), for I = 0, .., N — 1 . Up to a normalization 
factor, the kernel (4.28) is the Schrodinger functional we have been looking for. More 
precisely, we define 

Z = det(B N B )- 1 / 2 ( V +(+\T N [U N ,U ]\r )o ( ), (4.30) 
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thus discarding the determinants of the gauge field at the boundary. We are free to do so, 
because the boundary configurations are not integrated over. After all, the determinants 
have been introduced in order to cancel the Jacobians when the formula for the partition 
function is derived [4]. 

One may consider the Schrodinger functional as a functional of the fields at the 
boundary. In view of the continuum limit, we will now use the rescaled fields t\) c = 
(2/C) 1 / 2 ^ (cf. sect.l) with boundary values denoted as follows 

p +M x )\ Xo =o = />+(*)» p -M x )\ Xo =T = />'-( x )> 

(4-31) 

M x ) P -\x =0 = P-( K )> M x ) P +\ Xo = T = P'+( X ) ■ 

The generalization of the pure gauge Schrodinger functional (3.15) has the path integral 
representation 

Z[p' + , p'_, C ; p-,p + , C] = J B[t/,]D$]D[U] e~ S , (4.32) 

where the integration is over all quark fields at euclidean times xq with < xq < T and 
over the euclidean gauge field as specified at the end of section 3. The gauge part S g of 
the action S = S g + Sf has also been discussed there and the fermionic part Sf is given 

by 

3 

S f = a 3 V ip c (x)P_ f V a lk V k ip c (x) - U(x, 0)ip c (x + a 6) } 
xer fe=i 

T-a ^ 

+ ai J2 ^2M x ){^2^^ + m+-aA^fl} c (x) (4.33) 

xeTx° = a H 
3 

+ a z ^^ c {x)P + \^a lk V k ^ c {x)-U{x-aQ,Q)- 1 ^ c {x-aQ)^ _ . 
xer fe=i X °~ T 

Apart from the boundary terms, this is the Wilson action (2.6), again expressed in 
continuum notation. It should be emphasized that the whole boundary dependence is 
through the field components which fulfill the inhomogeneous Dirichlet conditions (4.31), 
the complementary components being absent. This can be understood from our derivation, 
since any Grassmann integral kernel only depends of tj and £ at the former and tj + and 
£ + at the later euclidean time. It can, of course, also be verified directly by inspection of 
the action (4.33). While this is easy to see as long as the lattice is present, it would be 
more difficult to establish in the continuum, because the difference between the fields at 
and near the boundary vanishes. 
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5. Continuum limit of the action 



As a first step in the investigation of the Schrodinger functional (4.32), one would like 
to know, whether it leads to a sensible classical theory in the continuum. In this section we 
thus consider the corresponding lattice action in the limit of vanishing lattice spacing a. 
To this end a gauge potential is introduced by parameterizing the link variables as 
follows 

U(x, fl ) = e a M x ) . (5.1) 

The classical continuum action S c is then obtained as the first term in the small a- 
expansion 

S[U,^,^] = S c [A,^^] + 0(a). (5.2) 
For the fermionic part of the continuum action, one finds 

S c f [A,j>,il>]= J dx J d 3 x^(x){j^D^ + m^(x) 



/Vx [^(x)P^(x)} - /Vx [^{x)P + ^{ 
Jo Jo 



(5.3) 



x = l 



Here, the covariant derivative is defined by = + A^, and the projectors P± = 
|(1 ± 70) are the same as in the preceding sections. The quark fields are ^-periodic in 
space and subject to the boundary conditions 

P+TP(X)\ X0 = = P+tt, P-^{X)\ X0 = T = p'_(x), 

(5.4) 

H x ) p -\ Xo =o = P-W> H x ) p +\ Xo =T = p'+( x ) ■ 

At first sight, the presence of the boundary terms in the action (5.3) is surprising. One 
might have expected to obtain the usual euclidean action, with the boundary condi- 
tions (5.4) imposed on the fields. At least, this is what happens in the pure gauge theory 
and the reader might be worried whether the boundary terms are a remnant of the lattice 
formulation. Stated differently, the question arises, whether it is possible to understand 
the boundary terms without reference to the lattice. 

To answer the question, one must first make precise what exactly one means by a 
sensible classical continuum theory. It will be assumed here, that the classical action is 
a functional on the space of smooth C°°-functions, which satisfy the boundary condi- 
tions (5.4). One may then set up a variational principle and look for the stationary points 
of the action. We will require the existence of smooth stationary points, which means that 
the corresponding solutions are required to be again in the above space of functions. 

Using these assumptions, one may investigate a fermionic continuum action with a 
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general ansatz for the boundary terms, 

S f [A,t(>,il>]= J dx J d 3 x^(a;)|7 Ai D Ai + m| 



+ / L d 3 x $(z)Ro(x)tl>(z)] + / L d 3 x $(x)R T (xW(x) 
Jo Jo 



(5.5) 



x = T 



The matrices i?o( x ) an d ^t( x ) parameterize the boundary terms and are left arbitrary 
for the moment. 

The variations SS and SS, defined through 

«g W = ] lm ^-» + "-d-^-*a, 

= , (5 ' 6) 

are distributions, evaluated on smooth test functions v and v. The test functions have 
to respect the boundary conditions on t\) and i.e. they are L-periodic in the spatial 
directions and satisfy homogenous boundary conditions 

P+v(x)\ Xo= o = P-<*)\ X0 =T = 0> 

(5.7) 

«(^-Lo=0 = ' *(*)P+\ X0 = T = ■ 

One then finds 

8S(v) = J dx J d 3 Kv(x)^ lt ,D t ,+m+S(x )P + R (K)+S(x -T)P_R T (K)^(x), (5.8) 
and, after integration by parts 

(•T pL 



SS(v) = J dx J d 3 x^(a;)|7 Ai [-V Ai + ^(aj)] +m 



(5.9) 

+ 6(x )(R (x) - jo)P- + S(x - T)(R T (x) + jo)P+}v(x) . 

According to the assumptions stated above, one requires the existence of smooth solutions 
tpd and tp c i of the equations 

SS(v) = 0, SS(v) = 0. (5.10) 

From eq.(5.8), one infers that tp c i must be a solution of (j^D^ + m)tp c i{x) = for all 
points x with < xq < T. The assumption of smoothness then implies that this must 
also be the case at the points xq = and xq = T, i.e. one must have 

P+M*)1>d{*)\ Xo =o = 0, P-R t (x)M*)\ Xo =t = . (5.11) 

14 



The same arguments, applied to tp c i lead to 

$ cl (x)(R (x)- l0 )P_\ X0=0 = 0, $ cl (x)(R T (x) + l0 )P + \ X0=T = 0. (5.12) 

Since one does not want to impose further conditions on the fields at the boundary, these 
equations are really conditions on the matrices i?o( x ) an d Rt(x), which must satisfy 

P+E (x) =0, (Eo(x) - 70 )P_ = 0, 

(5.13) 

P_E T (x) =0, (E T (x) + j )P+ = 0. 

The solutions are 

Eo(x) = -P- + P_r (x)P+ , R T (x) = -P+ + P+r T (x)P_ . (5.14) 

Here, To(x) and IV (x) are again arbitrary matrices. However, to give a non- vanishing 
contribution to the action (5.5), they have to anti-commute with 70. Such matrices are 
ruled out if one requires the action to be invariant under parity. 

In conclusion, one may say that the general ansatz (5.5) coincides with the fermionic 
action (5.3) as derived from the Schrodinger functional, provided the existence of smooth 
stationary points and parity invariance of the action are required. This is nice, because 
no reference to the lattice is needed, once the boundary conditions (5.4) are known. 



6. The free Dirac operator 

In the preceding section it has been shown how the action (5.3) can be understood 
without reference to the lattice formulation. An essential ingredient was the existence of 
a smooth stationary point, which implies that the decomposition 

S C [A, $ cl + v, i> cl + v] = S C [A, $ cl , i/> cl ] + S C [A, v, v], (6.1) 

does not contain linear terms in the fluctuation fields v and v. To get acquainted with the 
boundary conditions (5.7), it is helpful to have a closer look at the mathematical structure 
involved. Starting from the free action Sj[v, v], we first define the free Dirac operator. 
Motivated by potential applications we then establish its essential self-adjointness and 
determine its spectrum. 
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6.1 Definition of the free Dirac operator 

As usual, one would like to write the action S^[u,u] as a quadratic form in a suitable 
linear space of functions. At this point we observe that the boundary conditions (5.7) on 
v and v allow for the definition of two different linear spaces of smooth, spatially periodic 
functions u(x), namely 

H ± = {u\ P±u(x)\ Xq=0 = 0, P t u(x)\ Xq=t = 0}. (6.2) 
Equipped with the inner product 

def f T f L 4 
(u 1 ,u 2 ) n± = dx d 3 x Vti^ajJa^Wc, (6.3) 

Jo Jo a=1 

"H± are pre-Hilbert spaces. A smooth spinor field v is an element of "H + , while v defines 
an element G The free fermionic action can be written as a quadratic form in 

either "H + or "H-, 

S c f [v, v] = (v\ Dv) n _ = (£)t«t, v ) n+ , (6.4) 

where it is assumed that the operator D = 7^$^ + m and its formal adjoint act as 
follows 

D: H+^H-, : H-^H+. (6.5) 

Since both operators are elliptic (their leading symbols iip^j^ are invertible for / 0), 
this defines a structure, which is known in the mathematical literature as a two-term 
elliptic complex (cf. ref.[8]). 

We note that the eigenvalue problem for the operator D is not well-defined. However, 
according to Gilkey [9], one can obtain a first order elliptic boundary value problem if one 
defines an operator V as follows, 

V: U + ®U-^U + (BU-, ^= f (^ (6-6) 

The operator V thus acts on two-component vectors with the upper and lower components 
being elements of "H + and respectively. It is an elliptic operator which is also hermitian 
with respect to the induced inner product. We will refer to it as the Dirac operator in the 
following. 

A further notational item are the projectors 

which allow to characterize the linear space "H + as follows 

U = U + ®U. = {w\V + w(x)\ Xo=o = 0, V-w(x)\ Xo=T = 0}. (6.8) 
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The requirement that Vw be an element of "H entails the boundary conditions for the 
complementary components 

(0 O " rn)V-w(x)\ Xo=0 = 0, (d + m)V+w(x)\ Xo=T = 0. (6.9) 

With these additional boundary conditions, the Dirac operator squared, 

*=(7 A). <«•> 

also defines an elliptic boundary value problem. As will be shown below, the Dirac operator 
V and its square have unique self-adjoint extensions with purely discrete spectrum and a 
complete set of smooth eigenfunctions. 

6.2 Relevant questions 

If one wants to carry out a perturbative expansion of the Schrodinger functional, it is 
important to know whether the free quark propagator, i.e. the operator D~ l is well- 
defined. As will be shown below, the operator V 2 is essentially self-adjoint and strictly 
positive for any value of the mass. Eq.(6.10) then implies that this is also true for DD^ 
and D. It is therefore possible to define the propagator as follows 

D' 1 = D^DD*)- 1 . (6.11) 



In case dimensional regularization is employed, it is important to note that this statement 
remains valid for the (/-dimensional theory, provided translation invariance is not violated 
in the additional dimensions. 

Ultimately we will be interested in numerical simulations on the lattice. Current 
fermion algorithms involve the inversion of the matrix ^L at ^L at > i- e - the lattice version 
of D. The convergence rate of the matrix inversion is then determined by the so- 
called condition number, which is the ratio of the largest to the smallest eigenvalue [10]. 
Therefore, one has to make sure that the smallest eigenvalue is not too close to zero. For 
simulations in large volumes, this means that it is difficult to reach very light quark masses. 
The situation is different in small volumes, where the influence of the boundary conditions 
becomes important. In this case, one could envisage the simulation of light or even massless 
quarks, provided a smallest eigenvalue is induced by the boundary conditions. As will be 
explained in section 7, such simulations are required for the determination of the strong 
coupling at high energies, using the strategy of refs.[l,2]. 

We are thus led to investigate the continuum operator V 2 (6.10) and its spectrum in 
some detail. Close to the continuum limit, one expects that its lowest eigenvalues are well 
approximated on the lattice. It may be anticipated here that the boundary conditions (5.7) 
show their quality in inducing a minimal eigenvalue, given by (n/2T) 2 in the massless case. 
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6.3 Self-adjointness and spectrum ofV 2 

The operator V 2 is defined on functions w £ "H which satisfy the additional boundary 
condition (6.9). A considerable simplification is achieved by noting that, owing to spatial 
periodicity, V 2 is invariant under spatial translations. The eigenfunctions of this transfor- 
mation group are the plane waves e ipx , where the components pk of the spatial momentum 
p assume the values 2nnk/L, with integer numbers n^. The corresponding change of basis 
is given by the Fourier transformation 

w ( x ) = Ts y> ipx ™ P M, w p (x )= I d^e-^wix). (6.12) 
L P Jo 

In the following we restrict attention to the subspaces of definite momentum p, where the 
spatial derivative dk amounts to a multiplication by ipk- 

The operator V 2 is diagonal in Dirac space, so we may consider its action on each 
spinor component separately. A spinor component w p (xq)i is an ordinary function of xq 
which will be denoted by s(xq) in the following. One is thus led to study the differential 
operator — 8q + p 2 + m 2 on smooth functions s(x$) which satisfy either 

s(0) = 0, s'(T) + ms(T) = 0, or s(T) = 0, s'(0) - ras(0) = 0, (6.13) 

depending on whether s(xq) belongs to the subspace defined by the projector V+ (plus- 
components) or by V- (minus-components). From the literature (see e.g. [11], p. 237), one 
infers that this operator is essentially self-adjoint in either case. Moreover, its spectrum is 
purely discrete, all eigenvalues are non-degenerate and one may find an orthonormal basis 
consisting of its eigenfunctions. 

To find the eigenvalues and eigenfunctions, we write the corresponding eigenvalue 
equation 

(-do +P 2 +m 2 )s = X 2 s (6.14) 

in the form 

s"+p 2 s = 0, p 2 =X 2 -p 2 -m 2 , (6.15) 

which is simply the differential equation of a harmonic oscillator with frequency po. Taking 
into account the boundary conditions, the eigenfunctions are found to be of the form 
sinpo^o f° r the plus-components and sinpo( a; o — T) for the minus-components. In either 
case, one determines the allowed values for po as the solutions of the equation 

tanp T=-— . (6.16) 
m 

The case po = is special. The boundary conditions exclude a non-vanishing solution of 
s" = 0, except for the particular value of the mass m = — 1/T. In this case, there are 
linear eigenfunctions, which are proportional to xq (plus-components) or xq — T (minus- 
components). 
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Depending on the value of m, one has to distinguish 3 different cases for the solutions 
of eq.(6.16). For non-negative masses, all solutions are real. Since the eigenfunctions are 
odd functions of po one may restrict attention to the positive solutions. They lie in the 
intervals 

p (n )e K 2 ' / °Z M , no =0,1,2,... (6.17) 



For vanishing mass, the solutions are then given as the left limits of the intervals in (6.17), 
with n/2T as the lowest possible value. 

The second case corresponds to masses in the range > m > — 1/T. The solutions 
are again real and labelled as follows 

Po (n ) e T 2 ' ) , n = 0, 1, 2, . . . (6.18) 

In the last case, the mass is negative and below —1/T. Equation (6.16) admits real solu- 
tions in the same intervals as in (6.18) for no > 0. In addition, one finds a purely imaginary 
solution po(0) = iq, q > 0. One may continue to write e.g. sinpo^o f° r the correspond- 
ing eigenfunction, but one should be aware that this corresponds to the hyperbolic sine 
function isinhqxo. 

With these conventions, the complete set of orthonormal eigenfunctions of V 2 is given 
by {^n,z( a; )}> with 

/N f y/3(LT)- 3 f 2 ^[P + xo+P-(xo-T)]ei, if n = 0, m = -1/T, 

fn,l\ X ) = ) 

{ N(p ) e ipx [V+ sinpo^o + V- s'mp (x - T)] e/, otherwise, 



-1/2 



(6.19) 

Here, e\ , I = 1, .., 8 is a canonical basis in C 8 . The components n M run over all integer 
numbers if /x = 1, 2, 3, labelling the spatial momenta p = 27rn/Z/, and over all non-negative 
integer numbers if /x = 0, labelling the zero component po according to eqs.(6.17) or (6.18). 
The corresponding eigenvalues are given by 

X 2 n =p 2 + P 2 + m 2 , (6.20) 

each of them defining an 8-dimensional eigenspace. As suggested by the notation, the 
eigenvalues are all non-negative. The condition for the appearance of a zero mode may be 
derived from eq.(6.16), by substituting for po (6.15) and setting A to zero 



tanhyp^T^T=-^± m 



m 



(6.21) 
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The l.h.s. of this equation assumes values between and 1, while the r.h.s. is below — 1 
for positive masses and above 1 for negative values of m. A zero mode is thus obtained 
only in the limit of infinite negative mass. 

We may summarize the result obtained so far by saying that the square of the Dirac 
operator is an essentially self-adjoint, positive operator for all values of the mass. Its 
spectrum is purely discrete, with the lowest eigenvalue given by (n/2T) 2 in the massless 
case. Furthermore, its eigenfunctions (6.19) are smooth and define an orthonormal basis 
in the Hilbert space (i.e. the completion ofW). 

Concerning the Dirac operator itself, we note that its action on the eigenfunctions of 
V 2 reduces to the action of a hermitian matrix in each subspace of fixed n. One has 

8 8 

T> E W ^ n M = X) W 'l<Pn,l( X )> (6-22) 

1=1 1=1 

where w 1 £ C 8 is given by w 1 = Aw, with the matrix A defined as follows 

A= {a f)' A = W fe + (-l) no+1 VPo + ™ 2 - (6.23) 

This is a hermitian matrix with real eigenvalues, which are determined as the roots of 
A* (6.20), 

An,± = ±y / Po + P 2 +™ 2 - (6-24) 

Furthermore, since the trace of A vanishes, each eigenvalue is fourfold degenerate. 

Except for positivity, the operator V has the same properties as its square. In partic- 
ular, the set of functions (6.19) defines a core [12] which determines a unique self-adjoint 
extension of the Dirac operator. 



7. Concluding remarks 

So far, we have given a definition of the Schrodinger functional in lattice QCD. The 
formal continuum limit of the corresponding lattice action led to an acceptable classical 
continuum action, with a well-defined free Dirac operator. In this last section, we briefly 
discuss the renormalizability of the Schrodinger functional and comment on its application 
to the determination of the strong coupling constant. 
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7.1 Renormalizability of the Schr '6 'dinger functional 

The Schrodinger functional naturally leads to a field theory with prescribed boundary 
values for some of the field components. In general, quantum field theories with boundaries 
are afflicted by additional divergences. It is then a priori not clear whether the theory 
stays renormalizable. 

In his work on the Schrodinger representation in quantum field theory, Symanzik 
considered the case of massless scalar </> 4 field theory [13] (for an introduction to this 
paper see also ref.[14]). He found that the Schrodinger functional could be made finite by 
adding two new counterterms to the action. These are proportional to the local composite 
fields (j)do(j) and (f> 2 , integrated over the hyper-planes at xq = and xq = T. 

Following Symanzik, it is plausible that this behaviour is generic for a renormalizable 
quantum field theory. We thus have to look for possible boundary counterterms in the 
case of QCD. If we require invariance under parity, there are only two gauge invariant 
local composite fields of dimension 3 or less, namely 

■tptp, and tpjotp- (7-1) 

We therefore expect that the QCD Schrodinger functional is renormalizable after inclusion 
of at most two boundary counterterms, which are linear combinations of the composite 
fields in eq.(7.1). To confirm this expectation, it would be desirable to carry through a 
one-loop calculation with the QCD Schrodinger functional. This is currently undertaken 
in the framework of dimensional regularization, with the action (5.3) as starting point [15]. 

7.2 Application to the determination of a s 

It is our aim to apply the formalism to the determination of the strong coupling at high 
energies, following the strategy of refs.[l,2]. For the gauge fields, one could take over the 
formalism as developed for the the pure SU(3) gauge theory. In particular one could choose 
the same boundary conditions and the same definition of the running coupling. Concerning 
the quark fields, one would be tempted to set the boundary fields (4.31) to zero. If the 
renormalizability of the Schrodinger functional works out as expected, this implies the 
vanishing of the boundary counterterms [cf. eq.(7.1)]. The effect of the boundary is then 
of the order a and will thus modify the cutoff dependence of the theory. 

The Schrodinger functional serves to define a coupling which is assumed to run with L, 
the linear extent of the volume (one sets T = L). In order to trace its evolution from low 
to high energies, one has to cover a corresponding range of volumes. 

As has been shown in the preceding section, the boundary conditions force a minimal 
frequency into the system. In small volumes, or, equivalently, at high energies one therefore 
has a chance to simulate QCD with light quarks. As a first step, one could then try to 
simulate QCD with only massless quarks, thus staying with the bare coupling as a single 
free parameter. Of course, this is not quite true on the lattice with Wilson fermions, 
where chiral symmetry is broken explicitly. In principle, one must carefully tune the bare 
quark masses in order to find the chiral limit. However, it is plausible that the coupling at 
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high energies is not very sensitive to the precise values of the light quark masses, so that 
effectively only 1 parameter, the bare gauge coupling has to be tuned. 

It is difficult to overestimate the importance of a minimal frequency induced by the 
boundary. Otherwise, one would have to keep the quark masses at relatively large values 
so that a clean determination of the strong coupling might not be feasible. 

In sect. 6, only the free theory has been discussed. Of course, the gauge interactions 
will modify the spectrum of the Dirac operator and one expects the appearance of zero 
modes, corresponding to some particular gauge configurations. It will then be necessary to 
choose the background gauge field such that these gauge configurations become extremely 
unlikely. This will give an additional criterion for a good choice of boundary conditions 
for the gauge fields. At present it is not clear whether this criterion is satisfied by the 
background field which has been used for the simulations of the SU(3) pure gauge theory. 

I would like to thank M. Liischer for advice and many helpful discussions. 
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